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Abstract 

Within both slightly non-extensive statistics and related numerical model, a pic- 
ture is elaborated to treat self-similar time series as a thermodynamic system. 
Thermodynamic-type characteristics relevant to temperature, pressure, entropy, in- 
ternal and free energies are introduced and tested. Predictability conditions of time 
series analysis are discussed on the basis of Van der Waals model. Maximal magni- 
tude for time interval and minimal resolution scale of the value under consideration 
are found and analyzed in details. The statistics developed is shown to be governed 
by effective temperature being exponential measure of the fractal dimension of the 
time series. Testing of the analytical consideration is based on numerical scheme of 
non-extensive random walk. A statistical scheme is introduces to present numerical 
model as a grand canonical ensemble for which entropy and internal energy are cal- 
culated as functions of particle number. Effective temperature is found numerically 
to show that its value is reduced to averaged energy per one degree of freedom. 
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1 Introduction 



Time series analysis allows one to elaborate and verify macroscopic models of 
complex systems evolution on the basis of data analysis [1]. This analysis is 
known to be focused on numerical calculations of correlation sum for delay 
vectors that allow to find principle characteristics of the time series. Being 
traditionally a branch of the theory of statistics, time series analysis is based 
on the class of models of harmonic oscillator which are related to the simplest 
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case of the Gaussian random process [2]. But well known a real time series 
is relevant rather to the Levy stable processes, than the Gaussian ones being 
very special case [3]. Because the former are invariant with respect to dilata- 
tion transformation [4], the problem is reduced to consideration of self-similar 
stochastic processes. 

The simplest characteristic of a time series is known to be the Lyapunov 
exponent the largest of whose positive values yields predictability domain in 
the system behavior. A range of complexity in such behavior is determined by 
the Kolmogorov-Sinai entropy that equals to sum over positive magnitudes 
of whole set of the Lyapunov exponents 1 and can be reduced to the usual 
Shannon value in information theory. With passage to nonlinear system the 
probability p n of n-th scenario of the system behavior is transformed into 
power function p q n determined by index q < 1, so that the Kolmogorov-Sinai 
entropy should be replaced by the Renyi one 



K^-^—^Pl (1) 

H n 

Respectively, governing equation dp/de = —[3p q , (3 = const > describes 
probability variation with energy e = e n to derive the Tsallis distribution [5] 



p(e)oc[l-(l-g)/3e]^. (2) 

In the limit q — > 1, this distribution takes usual Boltzmannian form p(e) oc 
exp(— f3e) falling down exponentially fast in contrary to the power asymptotic 
of the generalized exponent (2). Physically, such a behavior is caused by self- 
similarity of non-extensive system [6] — [9]. 

As a result, the problem appears to study self-similar time series that present 
processes corresponded to power-law distribution type of Zipf-Mandelbrot 
ones [10]. This work is devoted to both analytical and numerical considera- 
tions of such type time series as a thermodynamic system. It is worthwhile to 
stress that our approach is principally different off the pseudo-thermodynamic 
formalism developed for consideration of multifractal-type objects (see [11], 
[12] and references therein). Indeed, if within latter formalism the role of state 
parameters play related multifractal indices, we introduce a set of effective 
state parameters, whose meaning is of true thermodynamic type (so, an effec- 
tive temperature is a measure of data scattering being related to the fractal 
dimension of the time series). 

The paper consists of two main Sections 2 and 3, the first of which is devoted 
to analytical consideration and the second one — to numerical study. We start 

1 We have in mind hyperbolic attractors. 
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analytical consideration with elaboration of a model which allows us to address 
a self-similar time series as slightly non-extensive thermodynamic system (see 
Subsection 2.1). Following Subsection 2.2 is based on calculations of both en- 
tropy and internal energy of the time series. As a result, t her mo dynamic-type 
characteristics such as temperature and entropy, volume and pressure, internal 
and free energies are introduced. Their testing for the model of ideal gas is 
shown to be basis for statistics of self-similar time series. Subsection 2.3 con- 
tains calculations of corrections to the ideal gas approach when external field 
and particle interaction are switched on. On the basis of Van der Waals model, 
we obtain the expressions for the specific heat and susceptibility as functions 
of the temperature. Subsection 2.4 is devoted to discussion of the physical 
meaning of the results obtained within the framework where the predictabil- 
ity of behavior of time series is mimicked as stability conditions of a non-ideal 
gas. We find maximal magnitudes for time interval and minimal resolution 
scale of the value under consideration. Subsection 2.5 shows that a tempera- 
ture governing time series statistics is exponential measure of a self-similarity 
index related to fractal dimension of the time series. 

Testing of the analytical consideration is based on numerical scheme of non- 
extensive random walk [13] whose stochastic equation and its solutions are 
treated in Subsection 3.1. As a result, we obtain non-trivial time series whose 
form is governed by a friction coefficient fixing related fractal dimension. In 
Subsection 3.2 we introduce a statistical scheme that allows us to consider 
modeled time series as a grand canonical ensemble for which we calculate 
entropy and internal energy as functions of particle number. In Subsection 3.3 
we find numerically effective temperature and show that its value is reduced 
to averaged kinetic energy per one particle of ideal gas. In concluding Section 
4 we confirm numerically exponential dependence of the effective temperature 
on the fractal dimension that was found analytically in Subsection 2.5. We 
discuss also multifractional and clustering properties of time series. 

Finally, several equalities needed in quoting are placed in Appendices A, B. 



2 Analytical study of self-similar time series 

The principle peculiarity of time series is that it evolves during large enough 
but finite time interval T < oo. On the other hand, the self-similarity means 
that a series relates to a fractal manifold characterized by the fractal dimension 
D. We show below that both properties pointed out are taken into account in 
natural way if we shall use non-extensive statistical mechanics where a simple 
combination of quantities T and D fixes the non-extensivity exponent q (see 
Eqs.(26), (31) below). 
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2.1 Description of time series as slightly non-extensive statistical system 



We start with consideration of ci-dimensional time series x(t;) related to the set 
{xj} of consequent values Xj = x(tj) of principle variable x(t) taken at discrete 
time instant tj = ir that we obtain as result of dividing a whole time series 
length T = Nqt by Nq equal intervals r. Following the ergodic hypothesis we 
shall imitate the time series x(tj) by the set of generalized coordinates {x n } 
supplemented by conjugated set {v„} of velocities that show jumping rates 
of coordinates Xj with the time variation. It is principally important to take 
into account that mapping of the time series x(tj) into the manifold {x n } 
is not mutually single valued because different terms x(tj) and x(t,) may be 
equal (similarly, it appears at mapping of velocity time series v(£j) into related 
manifold {v n }). This peculiarity is displayed first of all in that the number 
N of specimens of the time series is much more than the total number iV of 
terms of related manifold. In the following, just the latter number iV plays 
the role of the particle number of statistical system. Obviously, this system 
should be considered as grand canonical ensemble with variable number N n 
of particles in a state n. 

In the simplest case of Markovian consequence, one has the velocity v n = 
(x n — x n _i)/r. For more complicated series with m-step memory, the velocity 
magnitude is defined as follows: 



rn 



J2$?(m,n), 5i(m,n) = 



x (ra— m)+l 



— X 



(n-m)+(l-l) 



(3) 



Along the line of the ergodic hypothesis, the paradigm of our approach is to 
address the time series as a physical system defined by an effective Hamiltonian 
H = 7i{x„,v„} on whose basis statistical characteristics of this series could 
be found. If one proposes that series terms x n related to different n are not 
connected, the effective Hamiltonian is additive: 



N 

H = ^e n N n) e n = e(x n ,v n ). (4) 

n=l 

Physically, this means that the series under consideration is relevant to an 
ideal gas comprising of identical particles with energy e n and number N n in 
state n. Further, we suppose different terms of time series to be statistically 
identical, so that effective particle energy does not depend on coordinate x n : 
£(x n ,v n ) — > e(v n ). Moreover, since this energy does not vary with inversion 
of the coordinate jumps x n — x n _i, the function e(v n ) should be even. We use 
the simplest square form 
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Zn = ^Vl, (5) 

which is reduced to the usual kinetic energy for a particle with mass 1. With 
switching on an external force F = const, particle energy (5) becomes as 
follows: 



e n = \v 2 n - Fx n . (6) 

Finally, when time series has a microscopic memory, dimension of the delay 
vectors 5i(m,n) in the definition (3) needs taking m > 1. Moreover, if time 
series terms x m , x n with m ^ n are clustered, the Hamiltonian becomes 
relevant to a non-ideal gas with interaction w mn , m ^ n: 



1 N 1 

K = -Y.< + - 2 Y.^n. (7) 

n=l m^n 



To study a behavior of the time series as a whole one needs to fulfill summation 
over a set of states given by manifold {x„, v„} that is relevant to the system 
phase space. In so doing, it is convenient to pass to related integrations as 
following: 



V dx n dv ; 

{x„,v„} n=l ■ n=l 



e ^//n^=^- i n//dyndu n . ( 8) 



Here, the factorial takes into account statistical identity of the time series 
terms, A is effective Planck constant that determines a minimal volume of the 
phase space per a particle related to a term. The inverted factor 



/ V 2 \ - dN r ^9,v i -N 



eX 2d 

(9) 



_N(rA) d 
is caused by change of variables 



rescaled with respect to macroscopic length X being chosen to guarantee the 
conditions 

/dy n = l, |du n = l. (11) 
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According to Introduction self-similarity condition forces to use a statistics 
type of given by the Tsallis-type distribution (2). However, the latter is ap- 
peared to be inconvenient because the condition 



n 

takes place and definition of the internal energy reads: 



(12) 



n (1) £ 



(13) 



To take into account these constraints an escort distribution was proposed to 
use [14] 



V n = 



<1>, 



(14) 



In explicit form it reads as follows: 



V q {y n) u„} = < 



H{yn,u n }~E 



1-9 



at (1 - q) H{y $l- E < 1, 

(15) 



otherwise. 



Here, the partition function is defined by the condition 



N 



z = a' 1 n 



n=l 



l-(l-g) 



^{y», Un} - E 



1 

1-9 



dy n du n , 



(16) 



where < q < 1 is a parameter of non-extensivity, T s is energy scale. Internal 
energy E is determined by the equality 



A? 



n=l 



E = A/" -1 Yi // n{y m u n }T q {y n , u n }dy n du n , 



(17) 



and normalization parameter (1) is expressed by the partition function (16) 
in accordance with Eq.(88) in Appendix A. 

To check the statistical scheme proposed let us address firstly trivial case of 
time series x„ = const. Here, the particle energy e is a constant as well, so 
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that the Hamiltonian is 7i = Ne. The partition function Z = M^ 1 and the 
normalization parameter (1) = A/"^ 1-9 - 1 are given by inverted normalization 
factor (9), while the internal energy E = Ne is reduced to the Hamiltonian. 
Then, the entropy H = —a In A/", a = |(1 — q)dN ^ obtained according to 
definition (86) given in Appendix A is reduced to zero if only the normalization 
factor takes the value Af = 1. As a result, we find effective Planck constant: 



Our future consideration is stated on the assumption that the volume V = X d 
of rf-dimensional domain of the x n coordinate variation in dependence of the 
particle number N is governed by Levy-type law 



X d = x d N~. (19) 

Here, a; is a microscopic scale and a dynamic exponent z is reduced to the 
fractal dimensionality D of self-similar manifold [7] 



z — D. (20) 

Then, one obtains the following scaling relation for the phase space volume 
per a term of the time series: 



(21) 



In the case of Gaussian scattering, when D = 2, the minimal volume A d 
of the phase space does not depend on number iV of the time series terms. 
Such a condition approves our choice of the relation (19) for the whole volume 
V = X d as function of the number iV of time series terms. 



2.2 Non-extensive thermodynamics of time series as an ideal gas 



Calculations of main thermodynamic quantities of non-extensive ideal gas ar- 
rives at the following expressions for the partition function (16) and the in- 
ternal energy (17) (see [14] — [16]) 
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V N l(q) 
N\ 



dN 
2 



i + 



l + (l-q) — 



9) 



i 

1-9' 



dJV 
2 



(22) 



iV! 



f(l> t 



1-g 



dN 
2 



l + (l-<7) — 



1 + 



E 1 



1 , dJV 
l-q" 1 " 2 



(23) 



Here, <i-dimensional gas in volume V = X d is addressed and the notations are 
introduced 



2tiT s 
A 2 



7(9) 



2 J 



(24) 



to be determined by the Euler T-function. Combination of equalities (22), (23) 
with relation (88) yields explicit expression for the normalization parameter: 



(1) 



X dN i(q) 



6(l + a)i +1 
l-q 



1-9 
1-a 



(25) 



a = -(l-q)dN. 



In the limits 



1 - q < 2/d, AT > 1 



when 



(26) 



(27) 



7 (g) ^ [e(l-g)] 2 (1 + a) 
one obtains for the first of entropies (86): 



(28) 



H 



Na 



2(1 - a) 



In 



(29) 
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With accounting scaling relation (19), this expression takes the usual form 



H = N^nQ, G S (2xeT,)^g)" D (30) 
if the dynamic exponent is determined as 



z = D = ——. (31) 
1 — a 

Respectively, the internal energy (23) and the normalization parameter (25) 
read: 



„ dN/G\%^ .. (Gyr D-l 

e --{n) t - «.=y ' « s — ■ < 32 » 



The physical temperature is defined as follows [15] 



T = (l) q T s =(-Y T s (33) 

where the last equality takes into account the second of relations (32). This 
definition guarantees the equipartition law 



E = CT, C = cN, c = - (34) 



where the quantity 



is the specific heat. It is easily to convince that equations (30) — (33) arrive 
at standard thermodynamic relation 



OH 1 

be s r < 36 > 



Above used treatment is addressed to a fixed value of the internal energy E 
[16]. In alternative case when the principle state parameter is the temperature 
T, we should pass to the conjugate formalism [14]. Here, standard definition 
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F = E - TH 
of the free energy arrives at the dependence 



(37) 



F = -CT In 



T 

eT x 



(38) 



Then, the thermodynamic identity 



dF 



= -H 



dT 

yields the relation 



(39) 



H = Cln 



T 



(40) 



that plays a role of the heat equation of states. It arrives at the usual definition 
of the specific heat (cf. Eq. (35)) 



C = T 



dH 
~df' 



(41) 



Let us introduce now a specific entropy per unit time 



PjTT 

h=(dr)- 1 - m =r- 1 H 1 -r 
to be determined by a minimal entropy 



(42) 



H! = (D- 1) 
and a redundancy 



lis \/2nc'F — hi ( — 



(43) 



a i i , r \ D — 1 

r=-HeN), a=—. 

Dependencies on the scale x 



(44) 



h(x) = const — — In ^— ^ , r(x) = const 



(45) 
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notice that the system behaves in a stochastic manner [17]. 
Effective pressure is defined as 



9h , 
p=-r Tx (46) 

to measure specific entropy variation with respect to the time series scale. 
Then, we arrive at a mechanic-type equation of states 



px = D - 1 (47) 

being additional to the relation (40) of entropy to temperature. According to 
Eq. (47), definition of the pressure coefficient 



shows that it is fixed by the dynamic exponent (31) being the fractal dimen- 
sionality of self-similar time series: 



k — D — 1. (49) 

It is worthwhile to note that effective pressure (46) is introduced as derivative 
of the specific entropy h with respect to the microscopic scale x. That reduces 
a susceptibility of the time series to the pressure coefficient (48) inverted. 



2.3 Corrections to the ideal gas approach 



Let us focus now on effect of both external field and particle interaction whose 
switching on is expressed by equalities (6), (7). An external force F = const 
causes the second term in Hamiltonian (6) to arrive at the factor 



J ext 



sinh 



FX 
2T 



N 



(50) 



in partition function (22). Here, is a factor depended on the dimension d 
only and we put q — > 1 due to the conditions (27). According to the definition 
(86) the factor (50) yields the entropy addition 
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H. 



ext 



Na In 



sinh (ff) 



FX 
2T 



(51) 



where we suppress unessential term. With increasing homogeneous external 
field, the entropy H ext grows quadratically at FX <C T and linearly at FX ^> 
T. 

According to [18], cluster expansion of the particle interaction in Hamiltonian 
(7) results in additional factor in partition function (22): 



N 2 



Z int = 1 - — [ v + — 



w 



w 



S d J w(x)x d - 1 dx, S d 



2V V T 

2n d / 2 



(52) 



r(d/2)' v ~ d e 



where e is effective radius of the particle core. Relevant entropy addition 



Hi 



int 



aN 2 ( w 

v H 

2V V T 



(53) 



decreases monotonically with growth of the particle volume v. Increase of the 
temperature T causes the entropy decrease for attractive interaction w < 
and its increase in the case of repelling one (w > 0). 

Entropy additions (51), (53) arrive at total value of the specific heat (41) in 
the following form: 



C = -N — aN 
2 



\FX , 




— — coth 


v 2T ) X ] 


. 2T 





N w' 
2VT. 



(54) 



where formula (34) are taken into account. On the other hand, making use of 
equalities (42), (46), (51), (53) arrives at the pressure addition Sp determined 
by the equality 



5p ■ X = : 

a 



FX , (FX , 
coth f ) - 1 



N ( w 
+ a V\ V+ T 



2T \2T 

Then, with accounting (49) the pressure coefficient (48) takes the form 



(55) 



K 



a 

- + 1< 

a a 



1 + 



FX 
2T 



cosh (f^-) 



2^1 



(W 
V + T 



(56) 
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where effective density 



n = y = x~ d N a 



(57) 



is introduced. In the limiting case of low temperatures T <C FX/2, we obtain: 



C^^N-^(FX-nw), 



(58) 



a a 



k ~ 



/FX\* ( FX\ 
1 + ^— ) exp 



1 — a d 

Respectively, in opposite case T 3> FX/2, one has: 



+ (1 + d)an 



10 

T 



(59) 



C ~ -A^ + — 
2 2T 



nw — 



(fxt 

6T 



(60) 



1 + 



if 



FX 



4 V T 



+ (1 + d)an (v + 



(61) 



Thus, in the limits of both low and high temperatures, the influence of external 
field F reduces to hyperbolically decreasing addition (FX/T) n , n = 1, 2 to the 
specific heat C; accordingly, the pressure coefficient k gets the constant a/d 
only In a similar manner, the particle interaction affects hyperbolically on the 
specific heat, while the pressure coefficient takes the term being proportional 
to the particle volume v. 



2.4 Stability conditions 



It might seem the analysis of a time series would appeared to be the subject 
of more much study if both the interval r and the scale x take infinitely 
decreasing magnitudes (respectively, the number N tends to infinitely large 
values). However, we will show now that external influence and term clustering 
arrive at predictability boundaries in behavior of relevant system that are 
represented as instability of modeling thermodynamic system. Our analysis is 
stated on stability conditions 



C>0, k>0 



(62) 
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for the specific heat C and the pressure coefficient k given by equalities (54) 
and (56). 

In the simplest approach being modeled by the ideal gas in Subsection 2.2, we 
arrives at the natural restriction for the dynamic exponent: 



z = D > 1. (63) 

It means non-extensive dependence of the system volume (19) on the particle 
number N. On the other hand, in accordance with relation (31) the magnitude 
of the principle index (26) is limited by condition 



a < 1 (64) 

which restricts the rate of increasing effective density (57) with N. 

As is seen from expressions (54), (58), (60), the specific heat falls down mono- 
tonically with temperature decrease. To analyze such a behavior quantitatively 
let us consider the case of low temperatures T <C FX/2, when we can use the 
simplest dependence (58). Then, it is easily to see that at temperature less 
than critical magnitude 



T c = j [FxN^tt - wx~ d N a ) (65) 

the system becomes nonstable due to external force F and interparticle at- 
traction w < 0. If the number N is much less than critical value N c defined 

as 



Nc d ~^ a = ^x~ (1+d \ (66) 
F 

main contribution gets interparticle attraction w < 0, whereas in opposite case 
N 3> N c — external force F. With passage to more complicated case of high 
temperatures T 3> FX/2, when one follows to use temperature dependence 
(60), the physical situation does not change qualitatively. 

As show estimations (59), (61), the pressure coefficient becomes negative if 
microscopic scale x d is less than a critical magnitude 



^T^Pv"- <67 ' 

that is determined at attractive nature of interaction only (w < 0). If both 
parameters w and v takes non-zeroth magnitudes, the temperature values 
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are characterized by another minimal magnitude T v = \w\/v, lower which the 
microscopic scale x c becomes principle. This temperature is less than critical 
value (65) if elementary volume takes magnitudes lower than value 



v c = -n- 1 = -x d N~ a (68) 
a a 

where we take into account definition (57). 

In the limits v — > 0, F — > 0, boundary magnitudes (65), (67) arrive at appear- 
ance of the critical value 



a r = 



[l + d)(l + d 2 



Ad 2 



\ (i + d)(i + d 2 y 



(69) 



Relevant fractal dimension of the time series (31) grows monotonically with 
dimension of the time series d from the mag nitude D c = 1 + V2 ~ 2.414 at 
d = 1 to D c ~ d 2 at d ^> 1 (see Figure 1). As a result, the critical value 
(69) arrives at the upper boundary of the fractal dimensionality D c being 
more than topologic magnitude D t = 2d always (D c > D t ). Thus, we can 
conclude the upper boundary D c of application of our approach is appeared 
to be inessential at all. 





1 2 3 4 5 



Fig. 1. a. The dependence of the critical value (69) on the time series dimensionality, 
b. The same for related maximal value (31) of the fractal dimension of the time series 
(solid line) and topologic magnitude D t = 2d (lower line). The physical domain is 
shaded. 



2. 5 Effective temperature of time series 



We have considered above the simplest model which have allowed us to ex- 
amine analytically a self-similar time series in standard statistical manner. 
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Characteristic peculiarity of related equalities is a scale invariance with re- 
spect to variation of the non-extensivity parameter 1 — q which is contained 
everywhere through the parameter (26) that is related to the dynamic expo- 
nent z and the fractal dimension D of the time series according to Eq. (31). 
This invariance is clear to be caused by self-similarity of the system under 
consideration. For a given time series, the value D is fixed if it is addressed 
to a monofractal manifold and takes a closed set of magnitudes in the case of 
self-similar system relevant to a multifractal. 

In real time series the property of self-similarity can be broken, so that a 
dependence on the non-extensivity parameter itself could appear to be very 
weak. However, above introduced set of pseudo-thermodynamic characteris- 
tics of time series is kept as applicable and accustomed thermodynamic re- 
lations (34) - - (37), (39), (41), (46) and (48) can be applied to analysis of 
arbitrary time series. 

Main progress in our consideration is that time series statistics is governed 
completely by the temperature (33). With accounting the second of equalities 
(30) and rescaling the temperature unit T s into T sc = (27re) a T s we derive to 
the expression 



T 



r \ 



X) 



(70) 



Being independent of the number of terms N, the time series temperature 
shows exponential dependence on the index (26) located under critical mag- 
nitude (69). To establish a character of the power dependence on the ratio of 
the range X of the principle variable to the time interval r, it is naturally to 
choice measure units of the temperature in the following manner: 



T 

-*- fti 



A-a 



(27T) 



f) 2 



(71) 



Then, the expression (70) for the time series temperature takes the simplest 
form 



according to which the value T is the exponential measure of the fractal di- 
mensionality D of the self-similar time series. 

According to consideration given in Subsection 2.4, time series subjected to ex- 
ternal influence and term clustering is limited in predictability that is emerged 
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as instability of modeling thermodynamic system. On the basis of Van der 
Waals model, we have found minimal magnitudes (67), (68) for resolution 
scale x c of the value under consideration. Respectively, making use of the 
critical temperature (65) and the definition (72) shows that a time series is 
predictable if the microscopic time interval r takes magnitudes less than a 
critical magnitude r c . In the case N <C N c , where the boundary magnitude is 
determined by Eq. (66), we find 



(d A 2 D -i I d 2+d 1 
j\ w \) e^x + ^N"2Dd 2 . (73) 

In opposite case N 3> N c , one obtains 



_ l. 

(a „\ 2 d i i 
-F) e ~x2N^. (74) 

It is easily to convince that the entropy (40) takes positive values within the 
stability domain r < r c . 



3 Numerical study of time series 



In this Section, we are aimed to verify numerically the definitions based on 
expressions (3) — (5), (15) — (17), (26), (30) — (41). As above analysis has 
shown, one of peculiarities of the statistical system under consideration is that 
it is slightly non-extensive due to condition (27). Thus, we can put q — 1 in 
the following numerical consideration. 



3.1 Simulations of time series 



In our consideration, we shall follow the numerical scheme of non-extensive 
random walk [13] based on the discrete stochastic equation 



Xi+i = VrQ + (1 - 7t) + y/r£i 



(75) 



Here, discrete time ti = ir is fixed by integers i = 0, 1,...,N and minimal 
interval r; Q and & are additive and multiplicative stochastic sources normed 
with white-noise conditions (CiCj) = = friction coefficient 7 deter- 

mines parameter v = 7/(1 + 7) which fixes, in accordance with stationary 
distribution, non-extensivity parameter q = (2 — z/) -1 = (1 +7)/ (2 + 7). Mak- 
ing use of iteration procedure (75) arrives at stochastic time series, whose form 
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is shown in Figure 2 for total T = 500 and minimal r = 0.01 time intervals and 
different v values. As is seen from these series, at friction coefficient — 1 < 7, 
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Fig. 2. The form of modeled time series at T = 500, r = 0.01 and different v. The 
magnitudes of fractal dimension D are pointed out being determined in correspon- 
dence with J R/S'-method [20]. 

when the parameter v is negative, time series are relevant explicitly to the 
Levy flights. With passage to positive parameters 7, u, their growth arrives at 
gradual transformation of superdiffusion process into Brownian diffusion that 
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is related to the madnitudes 7 = 00, v — 1. As v increase, the relevant fractal 
dimension grows from the value D>lati/<0toD<2ati/<l. 

In Figure 3, we show velocity time dependencies v(t) obtained from origin time 
series x(t) according to definition (3) where the simplest non-clustering case 
m — 1 is taken. It is easily to see a velocity time dependence v(t) has much 
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Fig. 3. The form of velocity time dependencies v(t) related to time series x(t). 
Because R/ S-method gives overestimated values D R i S near upper boundary D = 2 
of fractal dimension, we apply also the values determined within box-counting 
method [20]. 

more rugged form in comparison with related time series x(t). As a result, 
fractal dimension for the velocity series v(t) is always more than its value 
for the initial series x(t) (see data in Figure 3). Moreover, it turned out that 
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R/ S-method and box-counting method (see [20]) give quite different values 
Dr/s and D^ c of fractal dimension. In our opinion, this difference which grows 
with passage from time series x(t) to velocity v(t) is caused by that related 
manifolds are more likely multifratals than fractals being corresponded with 
single value of the fractal dimension. 



3.2 Statistical consideration of simulated time series 

To study statistical properties of above time series we are stated on the er- 
godic hypothesis that supposes identity of averaging over both velocity time 
series v(ti), % — 1, 2, N and statistical ensemble {v n }, n — 1, 2, N, which 
describes velocity scattering in phase space. We determine such an ensemble 
dividing maximum interval of the velocity variation domain into N ^> 1 zones 
n— 1, 2, N, within of which velocities have mean value v n and small varia- 
tion Sv n . Then, the probability p n to hit the interval [v n — Sv n /2, v n + 8v n /2] 
and the relevant probability density function ir(v n ) are determined as follows: 

= (*>) 

where v n is a number of the time series specimens with mean value v n . Accord- 
ing to Figure 4 the probability density function 7r(t>) has the usual bell-shaped 




Fig. 4. The probability density function of the velocity distribution at 
different fractal dimensions (curves 1 5 correspond to parameters 

v = 0.85; 0.50; 0.35; 0.99; — 0.10; respectively, fractal dimensions are 
D = 1.81; 1.58; 1.45; 2.57; 1.45). 

form centered at the velocity v — 0. By this, non-monotonic variation of the 
tt(v) form in dependence of the fractal dimension D attracts attention to one- 
self: in particular, the most broad velocity scattering (curve 5) corresponds to 
negative value of the parameter v to display fat tails. 
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Because distribution (76) corresponds to grand canonical ensemble, related 
thermodynamic functions are determined by the n-state number 



N„ 



N 

N Pn = -u n 



(77) 



and energy (5). Thus, the total energy of ideal gas 2 reads: 

N N N 

71=1 



7V n=l 1 



(78) 



The plots of dependencies E(N) of energy (78) accompanied with related 
dependencies H(N) of entropy (87) are shown in Figures 5 for different values 
of fractal dimensions D. It is principally important, the energy E is directly 





Fig. 5. Dependencies of energy E and entropy H on the number N of effective 
particles (curves 1 — 4 correspond to fractal dimensions D = 1.50; 1.45; 1.58; 1.81). 

proportional to the particle number N, while the entropy H increases much 
more slowly with N. As for the velocity distribution tt(v), non-monotonic 
variations of both energy and entropy as functions of the fractal dimension 
take place, but main tendency is E decrease and H increase with D growth. 



3.3 Calculations of effective temperatures of time series 



The peculiarity of the self-similar system under consideration is that it is finite 
to be characterized with the following effective temperatures: 



6" 



dH 
~dE 



T- 1 = N 



D 



dH 
~dE 



(79) 



N 



Being determined at constant value of the fractal dimension D, the first of 
these is a function of the particle number N, whereas the second one depends 



Do not confuse with Hamiltonian (4). 
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on magnitude D to be determined at fixed value N. As show simple calcu- 
lations in Appendix B, above temperatures are connected by the following 
relation: 



T 



9 

N 



1 - 



1 



1 -q)D 2 \D - 1 



H D-l x ~ r 



D 2 



(80) 



Taking into account definitions (26), (31) and (34) accompanied with Eq.(40) 
where T s rescaled by factor N, one obtains 



T = 



e 

N 



Din 



T s J d D 



Making use of the data presented in Figure 5 gives for the first temperature 
(79) the dependence Q(N) shown in the left panel of the Figure 6. It is seen this 
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Fig. 6. Dependencies of effective and physical T temperatures on the num- 
ber N of effective particles (curves 1-4 correspond to fractal dimensions 
D = 1.50; 1.45; 1.58; 1.81). 

temperature takes very large values which increase directly proportionally with 
the N number growth. Such a behavior displays an intermediate character of 
the quantity G that plays rather a role of effective energy being complementary 
to the energy (78). On the other hand, usage of the relation (81) arrives at the 
dependencies of the temperature T(N) shown in the right panel of the Figure 
6. It is seen with increase of the fractal dimension to values D which are close 
to D = 2 the dependencies T(N) approach to constant magnitudes. These 
dimensions relates to large values of the parameter v < 1 which is connected 
with non-extensivity index q according to equality q = (2 — v)~ x . Thus, we 
find the temperature T takes values which are non-dependent on the particle 
number N in region q m 1 that is relevant to slightly non-extensive limit 
where above analytical consideration is valid. Due to such a behavior we can 
conclude the temperature T has usual physical sense. 

To confirm this conclusion we examine the equipartition law (34) rewritten 
with usage of the averaged kinetic energy (5): 
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(82) 



With this aim we compare the values (e) and T at different magnitudes of the 
fractal dimension D in Figure 7. It is seen the difference between above values 




Fig. 7. Dependencies of the averaged kinetic energy (e) and the physical temperature 
T on the fractal dimension D. 

is not exceed 10% if one does not take into account the points related to three 
smallest magnitudes of D where our analysis is not applicable. 



4 Discussion 



As shows above consideration, the quantity T defined with the second equa- 
tion (79) plays the role of the physical temperature that determines the energy 
(82) per one particle of slightly non-extensive ideal gas. According to connec- 
tion (72) this temperature should grow monotonically with increase of the 
fractal dimension, whereas Figure 7 shows non-monotonic dependence. As is 
seen from Figure 2, such a behavior is caused by strong falling down of the 
maximal variation X with D increase. To examine exponential relation in 
Eq.(72) we plot in Figure 8 the dependencies of the averaged special energy 
(e) and the temperature T on the fractal dimension D using the factor (X/t) 2 
as the scale of these quantities. It is seen with such a scaling above dependen- 
cies take approximately exponential form. Thus, main result of our analytical 
consideration is conformed numerically. 

It is worthwhile to pay attention to essential scattering of the numerical data 
presented in Figures 6 — 8. In our opinion, there are two main reasons of above 
scattering. Firstly, this is caused by multifractal nature of non-extensive time 
series generated according to recurrent procedure (75). To take into account 
this reason it needs naturally to consider a multifractional set whose fractal 
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<e>(X/z) 




2.4 D 



Fig. 8. Graphical check-up of the dependence (72): left panel — for mean value of 
kinetic energy (e); right panel — for the physical temperature T. 

dimensions are distributed as a function of the parameter q = (2 — v)~ l in 
accordance with the constraint 



D 



dN \-v 



-i 



v 



(83) 



that follows from Eqs.(26), (31). The second reason of the pointed out scat- 
tering is clustering of time series that is displayed essentially at small values 
of the parameter v (see Figure 2). This reason is taken into consideration by 
means of introducing effective interaction as it is done in Subsection 2.3. 



Appendix A: Main statements of the non-extensive statistics 



Our approach is stated on using entropy definitions 



H = a\n[exp q (H q )], H q = ln,[exp(if/a)]; a=-(l-q)dN (84) 

that are alternations of the usual functions logarithm ln(x) and exponential 
exp(x) with corresponding Tsallis generalizations [5] 



x 1 -? — 1 i 
ln q (x) = - _ , exp q (x) = [1 + (1 - q)x] . 

In explicit form, relations (84) appears as 



(85) 



H = a In Z, H q = 



1-5 ' 



(86) 
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Making use of the first equality (84) and the second formulae (85) and (86) 
shows that physically defined entropy H is extensive value being reduced to 
the Renyi definition (1): 



H = aK q =™\n(l) q . (87) 
Comparison of this result with the first of equalities (86) arrives at the relation 



(l) q = Z'^ (88) 
that ensures the normalization condition [5]. 



Appendix B: Relationship between effective temperatures corre- 
sponding to constant values of fractal dimension D and particle 
number N 



Within determinant representation, the first of definitions (79) reads as follows 
[19] 
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Making use of Eqs.(26), (31), (34) gives 
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- T. 

2 



(90) 



Respectively, rewriting relation (30) in the form 



S = N 



In(gN) + [D\ng + — \nN 
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(91) 
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one finds 
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A^ln N. 



Inserting Eqs.(79), (90), (92) in last relation (89), one obtains the result (80). 
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